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Abstract 

In order to investigate when neutrinos cease to oscillate in the framework of 
quantum field theory, we have reexamined the wave packet treatment of neu- 
trino oscillations by taking different sizes of the wave packets of the particles 
involved in the production and detection processes. The treatment is shown 
to be considerably simplified by using the Grimus-Stockinger theorem which 
enables us to carry out the integration over the momentum of the propagating 
neutrino. Our new results confirm the recent observation by Kiers, Nussinov 
and Weiss that a precise measurement of the energies of the particles involved 
in the detection process would increase the coherence length. We also present 
a precise definition of the coherence length beyond which neutrinos cease to 
oscillate. 

PACS number(s): 13.15.+g, 14.60.Pq 



1 



Typeset using REVTeX 



A rigorous treatment of neutrino oscillations requires the study of the processes in which 
neutrinos are produced and detected |]-f|. Following the first attempts to develop a 
proper use of the wave packet formalism, we have carried out detailed calculations in both 
quantum mechanics H and quantum field theory ||, with a quantitative derivation of the 
coherence length for neutrino oscillations. Due to the complexity of calculations, however, 
in || it was assumed that the sizes of the wave packets of the initial and final particles are 
all the same. 

Recently, Grimus and Stockinger || proved an elegant and very useful theorem that 
allows to simplify the wave packet treatment of neutrino oscillations. Taking advantage of 
this theorem, we have re-derived the neutrino oscillation formula in the general case in which 
all the particles involved in the production and detection processes have different wave packet 
sizes. Our new result confirms in the framework of a quantum field theoretical approach an 
interesting observation presented in |7j that an accurate measurement of the energies of the 
particles involved in the detection process leads to an increase of the coherence length for 
neutrino oscillations. 

Let us consider the flavor- changing process 



P I ^P F + C + v a 

v & + Di -»• Dp + 



(1) 



where Pi and Pp (Dj and Dp) are the initial and final production (detection) particles. The 
process ([]]) takes place through the intermediate propagation of a neutrino, which oscillates 
from flavor a to flavor j3 (here a,/3 = e,n,r). In the process (ffj), the production and 
detection interactions are localized at the coordinates (X P ,T P ) and (X D ,T D ), respectively. 

The form of the wave functions of the initial and final particles involved in the process 
(H|) is determined by how the initial particles are prepared and how the final particles are 
detected. In the following, we will assume, for simplicity, Gaussian wave functions, whose 
wave packet forms in momentum and coordinate space are given in || . The wave packets in 
momentum space are assumed to be sharply peaked around their average momenta, which 
are denoted by (f>k), where k = Pi, P F ,a, Dj, Dp, {3. The corresponding average energies 
(Ek) are given by (Ek) = \J (pk) 2 + m l, where is the mass of the k th particle. In order 
to make a realistic calculation, we will consider a different spatial width a x k for the wave 
packet of each particle involved in the process (|IJ). Let us emphasize that the localization of 
the particles does not require necessarily the action of a man-made apparatus, but can be 
determined by the environment in which the process (|l|) takes place. 

Following the method presented in ||, the amplitude of the process (fl|) can be written 
as (see Eq.(8) of 0) 



d^q 



qz _ m z _|_ %t 



Vp exp 



x J d 4 xi exp 
x / d 4 x 2 exp 



(2vr) 4 

-i (E P - q ) h + iipp- q) • xi 
-i (E D + g ) h + i (pd + q) • x 2 2 



-iq T + iq- L 



x\ — 2 vp • xi t\ + Tip t\ 



x% — 2v D -x 2 t 2 + Sd t\ 



(2) 



2 



where the index a denotes the neutrino field v a with mass m a , U is the neutrino mixing 
matrix, L = Xr> — Xp is the macroscopic distance vector from the neutrino production 
point to the neutrino detection point and T = Tp, — Tp is the macroscopic time interval T 
between neutrino production and detection. In Eq.(H) we have defined various quantities 
relative to the production process as follows (the corresponding quantities relative to the 
detection process are denoted with P replaced by D and a replaced by f3 in the subscript). 
Ep and ftp are given by 

E P = (Ep,) - (E Pf ) - (E a ) , (3) 
Pp = (m) ~ (Pp f ) - (Pa)- (4) 



The quantity a x p, given by 



111 

+ T^ + -^T, ( 5 ) 



a xP P^P! a xP F a xa 



can be considered as the effective size of the production process. Notice that the width of 
the wave packet of the most localized particle (i.e. the one with the smallest a x ) dominates 
in the determination of a x p. For example, in ]7| it is argued that the width a xe of the 
captured electron in the production process of 7 Be solar neutrinos, e~ + 7 Be — > 7 Li + v e , 
is much smaller than the widths of the wave packets of the other particles involved in the 
process. Therefore, in this case a x p ~ a xe (the formalism presented here can be obviously 
generalized for the case of more than one initial particle). 

In a quantum-mechanical framework (see ) the size of the production process 

determines the size of the wave packet of the propagating neutrino. On the other hand, in 
the approach adopted here the neutrino is treated as a "virtual" particle propagating between 
the production and detection vertices, whose properties are defined by the production and 
detection interactions in an equal way. Hence, any attempt to define a neutrino wave packet 
in the framework of the calculation presented here would lead to a definition of its size in 
equal terms of both a x p and <j x p>. Such an interpretation is unacceptable from a causal 
point of view. However, the present approach does not need an interpretation in terms of 
neutrino properties, because the neutrino is not directly observed, and we will show that, 
nevertheless, it leads to a rigorous derivation of the oscillation probability, including the 
effect of the coherence length. 

The quantities vp and Hp in Eq.(Q) are given by 

- _ 2 ( vpj v Pf v a \ 

Ef ^ p (^ + ^ + JL), (7) 

°xP 1 °xP F °xa 



where Vk = (f>k) / (Ek) is the group velocity of the k th particle. Finally, the definitions of 
U d and Vp are given by 

U D = J?{(PD F ), @D z ))M@f>)) ^ + 75) , (8) 

Vp EE (1 + l5 )Yv a ((P a )) J^((PP F ), iPPr)) , (9) 
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where Jf ((pp F ), (p*Pj)) and J?((pd f ), are the matrix elements of the weak currents 

of the production and detection particles. 

Carrying out the gaussian integrals over xi and x 2 , we obtain 



A af3 (L, T) oc J2 W ^(3a J -^y U D g2 _^ + ie Vp exp 



-iqoT + iq-L 



x exp 



ipp ~ q) 2 [{ e p ~ Qo) -(Pp-Q)- M 



An 2 



4a 2 P X P 



(Pd + q) 2 [(Ed + qo) ~ (Pd + q) • v D f 



4^ 2 d 



(10) 



where we have defined \p(p) = ^p(d) — v%m) anc ^ a pP{D) is defined by the relation 
cr a;P(D) cr P p(D) = 1/2. Notice that the integrations over x\ and x 2 in Eq.(fJ) did not pro- 
duce the usual (^-functions representing energy-momentum conservation at each interaction 
vertex. Indeed, in the wave packet treatment the energies and momenta of the particles 
involved in the process under consideration do not have a precise value, allowing for the un- 
certainty in energy-momentum conservation that is necessary for the occurrence of neutrino 
oscillations (see [^HMl ) • 

From Eqs.(||) and (|7|) one can see that < Xp(D) < 1- Hence, it is natural to ask what 
happens if Xp = and/or Xd = 0. Let us consider, for simplicity, the case of only Xp = 0. 
This situation could arise, for example, if only the initial particle Pi in the production 
process is observed with a good localization at Xp and is described by a wave packet, 
whereas all the particles in the final state are described by plane waves. In this case we 
have a xP = a x p n vp = vp I and Sp = vp , leading to = 0. From Eq.fl2|) one can see 
that Sp = Vp implies exact energy conservation in the production process (5(qo — Ep)). If 
Vp 7^ 0, there is also exact momentum conservation in the production process [5 3 (q — pp)) 
and neutrino oscillations do not occur because exact energy-momentum conservation in the 
production process can be satisfied only by one of the neutrino mass eigenstates, excluding 
the coherent production of more than one neutrino mass eigenstates that is necessary for 
neutrino oscillations. This fact is also clear in coordinate space: if vp ^ only the trajectory 
of Pj is known, which is not sufficient to determine the position Xp and the time Tp of the 
production process. Hence, the distance L and the time interval T are not defined and 
oscillations are not observable. On the other hand, in the special case Vp = there is 
no exact momentum conservation in the production process (although energy is exactly 
conserved) and neutrino oscillations are observable. Physically this situation corresponds to 
have the initial particle Pi at rest, in such a way that Xp is known and the distance L is 
defined, although Tp and T are not known (similar scenarios are discussed in P||). However, 
as stated above, the calculation presented here is based on the assumption that all the 
initial and final particles in the process (P are observed and their wave packets are sharply 
peaked around the corresponding average momentum. Therefore, the case Xp = (and/or 
Xd — 0) corresponds to a different physical process from the one under consideration, which 
can be discussed modifying the calculation presented here in an appropriate way. Under 
our assumptions, Xp and Xd are dimensionless quantities of order 1CT 1 . For example, if 
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the production process is pion decay at rest, we have Xp 

(ml - ml) 2 /(ml + m 2 ) 2 , i.e. ^~7x 10~ 2 for a = // and v a — ± iui u — c. 

Resuming the calculation of the amplitude fllQf) , we perform the integral over q using the 
Grimus-Stockinger theorem M: 



d 3 q 



iq-L 



<p(q) e iq - 
ql-cf + ie 



L— >oo 



2tt 2 



(p(q a i) e 



iq a L 



(11) 



with L = \L\, £ = L/L and q a = y q° 2 — m\. This theorem is valid for a function <f> which 
is differentiable at least three times such that itself and its first and second derivatives 
decrease at least as 1/q 2 as |g| — > oo. This is precisely our case. The remaining integration 
over q° can be done with a saddle-point approximation at q° = E a , with E a given by the 
relation 



E a 



Pa 



■VP 



p a - E a £■ V P ) (E P - E a ) - (p P - p a t) • V P 



pP 



o- 2 pXp 



(p a - E a £■ Vd) [(E D + E a ) - (p D + Pa t) 



pD 



o- 2 d Xd 







(12) 



where p a = q a (E a ) = sjE 2 — m\. The quantities E a and p a can be interpreted as the 
effective energy and momentum of the a th neutrino mass eigenstate propagating between 
the production and detection vertices. In order to simplify our discussion, in the following 
we omit explicit terms of the order m^/E 2 , which are negligible for relativistic neutrinos. 
Thus, for L — > oo we have 



A a p(L, T) OC -jr ^ Ka U da A A= exp 
a V a 

where v a = p a /E a , 

A a = U D (E al -Pat-llVp, 



-iE a T + ip a L - S a (E a 



1 

2^ 



+ 



n 2 a 2 
°pP °pD 



+ 



■ Vp 



°Ip^p 



+ 



■ v D 



S a (E a 



(/PP - Pat) (Ep - E a ) - (p P - Pat) • Vp 



+ 



Pd+P. 



at) 



An 2 

^°pD 



+ 



4a 2 P \ P 
(Ed + E a ) - (p D + p a t) ■ V D 



1 (L-VgTY 

2 v 2 a n a 



+ 



mi 
E 2 



(13) 

(14) 
(15) 



(16) 



The probability of the process ([]]) is proportional to \A a p(L, T)\ 2 , but in a practical 
experimental setting L is usually a fixed and known quantity, whereas T is not measured. 
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Therefore, the oscillation probability at a given distance L is given by the time average of 
\A a/ 3(L, T)| 2 , which leads to 

P aP (L) oc ±- 2 Y,AaAlW aa U Pa U ab U; h -/=J= exp [S a (E a ) - S b (E b )} 



x exp 
x exp 



, v a VL a + v b tt b , . 

(£ a - E b ) /0 ; , - (p a -pb))L 



V a V b (^a + ^6 



(v - v b y 



2 v y b (n a + n b ) 2 



- (E a — E b 



Q r ,Q 



(17) 



The overall factor 1/L 2 represents the geometrical decrease of the neutrino flux. Notice that 
the probability ( |P7|) depends not only on the modulus L of the distance, but also on its 
direction £, which is contained in A a , Q a and S a (E a ) (see Eqs.(fl4])-(|r6|)) and determines E a 
through Eq . (|12"D . This is due to the fact that the wave packets of the external particles are 
assumed to be known and the term exp [—S a (E a ) — S b (E b )} in Eq.flTTI) guarantees that the 
probability of the process ([]]) is not negligible only when both pp and pp> are aligned with I 
within the uncertainty allowed by the sizes of the wave packets. 

Since we are concerned with relativistic neutrinos, we approximate 



E ~ 

J -'a — 



(18) 



with E and p determined, respectively, by the expansion of Eq.(^) at zeroth and first order 
in powers of m 2 a / E 2 . Equation flOp leads to the approximations p a ~ E + (p — 1) m 2 ./2E, 



1 - m 2 a /2E 2 and Vt a 



2uxt 2 , with 



2 , 2 
<J xP i* °xD 5 



UJ = 1 + 



a 



1 - 



v P 



a 



xD 



1 - 



VD 



(J: 



A 



D 



(19) 
(20) 



From Eq.(|T6D one can see that in the relativistic approximation S a (E a ) is minimum for 
Pp ~ — po — Ei and Ep ~ —Ep, ~ E. These are the values of the kinematical parameters 
for which the probability ( |P7D of the process (||) is not negligibly small (without consid- 
ering the geometrical suppression factor 1/L 2 ). Hence, energy-momentum conservation is 
approximately satisfied in order to guarantee the observability of the process ([!]) and this 
approximate conservation determines the value of the effective neutrino energy E. From 
Eq.(|l6D one can also see that p is a dimensionless quantity of order unity. On the other 
hand, to could be rather large if \p(D) is small. For example, if the production process is 
pion decay at rest with a — p,, Xp> ~ 1 and o x ~ cr xP , we have to ~ 10. 

In the relativistic approximation one can factorize out of the sum over the mass eigenstate 
indices a and b in Eq.(|T7p all the quantities that do not vanish in the zeroth order of 
the expansion in powers of m 2 a / E 2 {e.g. A a and S a (E a )). Therefore, in the relativistic 
approximation for the flavor-changing probability we have 



6 



P a p(L) = ^Ul a U Pa U ab Up h exp 



a.b 



-27T2- 



L 

T OSC 



T coh 
^ab 



- 2tt 2 lu P 2 



T osc 
^ab 



(21) 



with the oscillation lengths L°l c and the coherence lengths L c ^, for a^b, given by 

4vr£ 



- coh 
J ab 



Am lb ' 

,— 2E 2 
2V2uj — — 5-7 a x 

\&Kb\ 



(22) 
(23) 



where Am^ = — m 2 . 

The transition probability (|21|) has the same form as that given in Eq.(23) of |J, which 
was obtained in the framework of quantum field theory with wave packets. However, the 
coefficients that appear in the expression of the coherence length and in front of the third 
term in the exponential in Eq. fl2i"|) are different from those of || . There are two reasons for 
these differences: 1) in we assumed the same spatial width for all the wave packets of the 
external particles involved in the process (|IJ), whereas here these widths can be different; 2) 
in || we integrated over the momenta of the final particles, whereas here we assume these 
momenta to be measured. 

The first term in the exponential in Eq.Q2~T|) is the usual oscillating phase which gives 
rise to neutrino oscillations. The second term causes a quadratical decrease of P a p(L) 
with the distance L and determines how far the oscillations take place. For L > this 



term suppresses the interference of the neutrino mass eigenstates v a and vt,, leading to the 
disappearance of the oscillations due to L° s f for L ^> L^ 1 . If L ^> for all aj^b, all 
the oscillating terms in the probability (pl|) are suppressed, leading to a constant transition 



probability P a/3 = J2 a l^aaf \M/3a\ 2 - The third term in the exponential in Eq.(^), which 
is due to the time integration, i.e. to the lack of time measurements, implies that the 
interference terms are also washed out if a x is larger than the neutrino oscillation length. 

Apart from the factor 2\/2u;, the coherence length given in Eq.(P3"D is similar to that 
obtained by physical intuition in |l|-|3| and in a quantum mechanical wave-packet treatment 
in 0] . But there is a very important difference between our result and the previous formulas 
for the coherence length. The coherence length is usually defined as the distance at which 
the separation between the wave packets of different massive neutrinos propagating with 
group velocities v a ~ 1 — m 2 a / E 2 is equal to a x , i.e. 



t coh 

^ab = 



Vb 



2E 2 
\ Am lb 



(24) 



where a x is the width of the propagating neutrino wave packet, which depends on the 
neutrino production mechanism. However, our calculations show that the proper definition 
of the a x that determines the coherence length must also include information on the neutrino 
detection mechanism. We presented a rigorous definition of a x in Eq.(|i~9D. Defining a p with 



the relation a x a p 



1/2, we have 



Hence, it is clear that the smaller of a p p and g p b dominates the value of a p . If a p p, is much 
smaller than a p p, then a p ~ o p p> and a x ~ u x d- Therefore, a precise measurement of the 
momentum of all the particles involved in the neutrino detection process^ implies a very 
small a p ~ <t p d, leading to a very large a x ~ a x p> and a very large coherence length. Thus, 
our formulation with wave packets confirms the observation in and provides a simple 
method to estimate the a x that determines the coherence length through Eq.(^3|): o x ~ 
V a xP + a "xDi where a x p and <j x d are the estimated sizes of the production and detection 
processes, respectively. 

The presence of the factor 2\/2uj in the expression of the coherence length can imply 
a non- negligible increase of L c ° b h with respect to the usual definition (|2l|). For example, as 
we have already mentioned, if the production process is pion decay at rest with a = fi and 
\d ~ 1, we have lo ~ 10, which gives 2y2u ~ 10. Since to is large if Xp(D) is small, i.e. 
if all the particles involved in the production (detection) process have a very low velocity, 
this increase of the coherence length is due to the fact that the overlap of the wave packets 
of slow particles last longer and there is more time available for the production (detection) 
process to emit (absorb) coherently the intermediate superposition of massive neutrinos. 
Indeed, one can estimate, for example, that the time Atp available for the production 
process is of the order of <J x p/\vp\ and, if \vp\ is small, \z> ~ 1 and a x ~ a x p, we have 
uj ~ l/|?Tp| 2 ~ (Atp/a x ) 2 . Hence, one can see that a small \vp\ implies a large Atp, a large 
uj and a large coherence length. 

On the other hand, it must be noticed that a large value of uj increases the contribution 
of the third term in the exponential in the transition probability (^Tj) and could enhance its 
suppression of the interference terms if a x ~ L°l c . This is due to the fact that when the 
time available for the production (detection) process is longer than L°\ c , the interference of 
the neutrino mass eigenstates v a and v h is washed out. For example, if uj ~ (Atp/a x ) 2 as in 
the example above, we have u {a x j L°l°) 2 ~ (Atp/L™ b c ) 2 and the corresponding oscillating 
term in ( |2TD is suppressed if Atp ^> L°^ c . 

In summary, the Grimus-Stockinger theorem considerably simplifies the wave packet 
treatment. The neutrino oscillations due to L°| c disappear when the distance L between 
neutrino production and detection is much larger than the coherence length L^ h given by 
Eq.(^3D with Eq . ([19|) , which is the most complete expression derived so far. The authors 
of |7j pointed out that one can observe neutrino oscillations beyond the usually defined 
coherence distance by a precise measurement of neutrino energy. Our calculation confirms 
their claim in the framework of a quantum field theoretical calculation. Another result of 
our calculation is that, because of the lack of time measurements, the neutrino oscillations 
due to L°| c can be observed only if a x is smaller than L^ b c . Otherwise, the interference of 
the mass eigenstates v a and i>b is washed out. 



1 It is necessary to measure accurately the momentum of all the particles because from Eq.([|) we 
have a 2 pD = a 2 pDj + a 2 pDp + a 2 pp . 
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